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Abstract 

We study the linear stability of a circular orbit in a two-electron atom, with 

the inclusion of retardation and self-interaction effects. We calculate all the 

eigenvalues of the linear stability of the circular orbit, expanded up to third 

order in (v/c). We discuss how the linear stability results can be used to 

construct a resonant normal form constant of motion, which is important for 

the short time stability of the orbits (e. g. emission of sharp spectral lines). 

We calculate the magnitudes related to some of those orbits and compare to 

the quantum atomic results. 
PACS numbers: 31.15.Ct,03.20.+i,05.45+b 
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I. INTRODUCTION 



We have recently proposed that the retardation of the Coulomb interaction could stabilize 
some circular orbits of the helium atom by the existence of an extra constant provided by a 
resonant normal form ffl. This extra constant would return the coulombian radial instability 
to the neighborhood of the original circular orbit in a short time scale. In this letter we 
study in full the linear stability of a circular orbit of a two-electron atom with the inclusion 
of retardation and self-interaction. The linear stability is one of the ingredients in the 
construction of the resonant normal form. Another work will be dedicated to the detailed 
construction of the resonant normal form, which we briefly discuss here in appendix C in an 
intuitive way for completeness. 

A circular orbit is one where the two electrons are in the same circular orbit and in phase 
opposition, that is, along a diameter ||. The center of the two-electron atom is supposed 
to have a positive charge of value Ze, (Z > 1/4) and we henceforth call it the nucleus. The 
linearized dynamics about a circular orbit has one unstable direction, one stable direction 
and ten neutrally stable directions 0. For an infinitely massive nucleus, a two-electron 
atom has a six-degree of freedom Hamiltonian system with only four independent constants 
of motion || . (Namely, these constants are the energy and the three components of the total 
angular momentum). Because there are only four constants, the Coulombian dynamics in the 
neighborhood of a generic circular orbit can in principle be unstable (by lack of constants), 
and actually is unstable in this case f2|. Here we show that an extra complex constant can 
exist in the neighborhood of a discrete set of circular orbits, when one includes retardation 
effects. This provides two extra constants, and together with the four global constants of 
the two-electron atom, this can stabilize the dynamics in the region inside the resonance 
islands, in a short time scale. This extra constant is provided by a resonant normal form 
(adelphic integral and requires a resonance to exist 

Historically, the understanding of the electrodynamics of a charged particle interact- 
ing with its own electromagnetic field PHI3| came very late. A classical solution to the 
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self-interaction of a charged particle of very small radius is described by the Lorentz-Dirac 
equation of motion, henceforth called LDE. The derivation of this equation was first pre- 



sented in a lecture by Lorentz in 1906, and first published in 1909 PJIIf . Lorentz's theory 
had many difficulties ||13|| , and a major progress came only in 1938, when Dirac |10[ produced 
a covariant derivation without mention to the structure of the particle. Dirac was also the 
first to recognize and understand the runaway solutions to the LDE |]13| . Dirac's theory still 
suffered from an arbitrary mass term ||11|| . In 1948, the regularization approach invented to 



treat the Lamb shift, offered a satisfactory solution to the divergent mass renormalization 



In this work we consider the electrodynamics of pointlike charges with a renormalized 
mass, as described by the Lorentz-Dirac theory ]ll]-0. 

Another late development of Maxwell's theory was the work of Page (1918) |14| on the 
expansion of the Lienard (1898) and Wiechert (1900) formula. This formula is complicated 
because of the retardation constraint, and one way to convert it into a useful differential 
equation is to develop the constraint in a Taylor series. This was done by Page up to the 
fifth order in 1918 [[TJJ], who also explored this formula in connection with self-interaction. 
We henceforth call the expansion of the Lienard- Wiechert interaction the "Page series". 
Truncated to second order in (f/c), the Page series describes a Lagrangian interaction. 



This Lagragian is the Darwin Lagrangian |13],[15],[T6fl , which introduces the first retardation 



correction to the Coulomb dynamics. The Darwin Lagrangian is used in quantum mechanics 
to produce the Breit operator, which is the generalization of Dirac's equation for two-electron 
atoms, correct to second order in (v/c) [17],[18[]. The Breit operator has been used successfully 



to describe the spin-orbit coupling and fine structure of helium [|18"l , |19f . The third-order term 
of the Page series is dissipative, and of course not Lagrangian. 

The most studied case of a two-electron atom is the helium atom. The Coulomb dynamics 
of helium is not completely integrable, because only three constants of motion in involution 
exist for the dynamics 0,0. This non-integrability of helium appeared historically as 
a hindrance for the early quantization attempts of the Coppenhagen school Quite 
recently, because of the renewed interest in periodic orbit quantization, much has been 



learned about the Coulomb dynamics in helium |2^j2^]. For example, the Coulombian 



radial instability is now known to lead to self-ionization, after a long-term chaotic transient, 
for most initial conditions ||23|| . In section IV 111 we discuss how the six constants can stabilize 



the dynamics in the small region inside the resonance islands, in a short time scale. 

The introduction of the first retardation correction to the Coulomb interaction modifies 
the dynamics in a quantitative and in a qualitative way, as far as this work is concerned: 
The Darwin Lagrangian is rotationally invariant, which generates an angular- momentum- 
like constant of motion according to Noether's theorem [POR . This constant is the angular 
momentum of the electrons, with relativistic correction, plus the angular momentum of the 
electromagnetic field. There is also an energy-like constant (because the Darwin Lagrangian 
is time independent). These are small perturbations of the four constants of the Coulomb 
dynamics, and this is the quantitative change. The qualitative change is the existence of 
the extra analytic constant in the neighborhood of some orbits, and we stress that this is a 
genuine nonlinear effect, because the frequencies of the linearized dynamics with retardation 
depend nonlinearly on the orbit's frequency . This extra constant appears only after one 
includes the retardation effects, which unfold a degeneracy of the Coulomb dynamics. In 
this sense, we have included the retardation because it is absolutely necessary, and not to 
attain some better precision. 

The main motivation of our stability studies was to investigate the possible existence 
of an extra analytic constant of the motion in the neighborhood of some special circular 
orbits, as described by the Darwin Lagrangian. According to standard normal form theory 
0-0], the existence of an extra analytic constant of the motion in the neighborhood of 
an equilibrium point requires that the eigenvalues of the linearized dynamics satisfy some 
resonance condition. Of course, an extra constant would require a new, special resonance 
to be satisfied. We use the results of to find where the extra analytic constant can exist 
and we show how to construct this extra constant in appendix C. 

In this paper we accomplish two things: First we perform the complete stability analysis 
of a generic circular orbit, including retardation and self- interaction up to third order in 



(v/c) , which is done here in full detail for future reference for a generic two-electron atom.. 
A lot of the work goes into solving for the linear stability problem, which we call "variational 
dynamics" or sometimes tangent dynamics and we develop a systematic method to handle 
the variational equations with inclusion of the retardation and self-interaction corrections. 
As an application of the linear problem that we just solved, we explore some of the very 
interesting resulting orbits about which an extra constant is possible. We show that a 
degeneracy determines which eigenvalue enters the resonance condition to produce orbits in 
the correct atomic magnitude. We discuss mainly the cases of helium and the Li + ion. In 
this work we do not prove which resonances do stabilize the orbit in a linewidth time scale, 
which is the subject of another work. Since this requires moving to the rotating frame and 



performing long normal form calculations, we chose to devote a separate work to it |p4 |. To 
make the present work intuitively complete, we included in appendix C the main ideas of 
how to construct the extra constant, that can possibly secure nonlinear stability. We also 
discuss (in Section II) how the resonant terms of the normal form can produce a deformation 
of the Coulombian orbit in agreement with a sharp line. 

Some cautionary notes are in place before one starts reading this paper: First, we are 
studying dynamics in the neighborhood of a circular orbit, which is also a periodic orbit. 
One should not think though, that we are doing periodic orbit quantization fHH. It is known 
that in general EBK can not be applied to generic orbits of Coulombian helium because of 
the unstable direction |J. Since we are showing that the inclusion of retardation stabilizes a 
discrete set of orbits, it would be of interest for future research to apply EBK quantization 
in the neighborhood of those stable orbits, which would be a nice new application of normal 
form theory. Of course one would have to start the quantization not with the Coulomb 
Hamiltonian but with the Hamiltonian resulting from the Darwin Lagrangian, which includes 
retardation effects. |25 



The dynamical system with retardation and self-interaction has nothing to do with the 
corresponding quantum system in principle (unless otherwise proved in the future). It 
is nevertheless very interesting to compare our stability results to the quantum results. 



This comparisom immediately exposes the relevance of the results to atomic physics and 
suggests what further results should be sought by use of nonlinear dynamics. The first 
striking coincidence is that our dynamical orbits have energies around the correct atomic 
magnitudes and the emited frequencies agree very well with some frequencies of the spectra 
of helium and Li + . The essential novelty in the work is that we are discussing for the first 
time a dynamical system that can emit a sharp line (to emit a sharp line the orbit must 
be stable in a time scale of some 10 6 turns, the type of stability that resonant normal form 
can provide). Of course this would suggest a deeper interface between quantum mechanics 
and classical electrodynamics with retardation, and we hope that this interface will become 
more transparent as we go along with the research. We do not know at present if it is 
possible to approximate all the lines in the spectra of helium and Li + in this way. In our 
approximation, we retained only the minimal amount of terms absolutely necessary to arrive 
at the results. It comes nevertheless as a surprise that nonlinear dynamics can produce such 
a sensible approximation to some of the sharp frequencies of the spectrum of helium and also 
to stabilize circular orbits in the correct atomic energies. Last, in this work we consider only 
the interaction of pointlike spinless particles, and we are not including spin in the dynamical 
system. 

This paper is organized as follows: In part || we discuss the electromagnetic formulas. 



In part pj we consider the Coulomb stability problem, which is order zero of the Page series 
formula. In part |V| we discuss the inclusion of the second-order terms and consider the 
resonance condition necessary for an extra constant and compare some results. In part V 
we consider the influence of the radiative terms, in part VI we compare our results to the 
atomic results and in part [VII| we put the discussions. 



II. ELECTROMAGNETIC FORMULAS 

In this work we include the self-interaction effects as described by the relativistic Lorentz- 
Dirac equation (LDE) with a renormalized mass JT0|-|r5! . The LDE equation for an electron 
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of charge — e can be written in the convenient, noncovariant form as 



j t ( 1 M e ± e ) = T + F exU (1) 



where x e is the electron velocity, M e is the renormalized electronic mass [11] and 7 = 



1/y 1 — (|x e |/c) 2 . In ([!]), F ext is the external force acting on the electron and T is the 
radiation reaction force. For circular orbits, the lowest order term of T in powers of (v / c) is 
T = x e . The next correction to T is of order (v/c) 5 for a circular orbit |Tl . 



We now introduce, for later use, the expansion of the retardation constraint of the 
Lienard-Wiechert interaction, (the Page series [[14|]). Let x be the position of a charge 
q, and (3 its velocity vector divided by c. The formula for the electric field caused by this 
charge q at a point p is 

n , „ r (|/3| 2 -3(n-/3) 2 )n $ (n-flfi , 
p V +gl 2r 2 2rc 2rc 1 

where n is the unit vector pointing from the charge to the point p, r is the distance from the 
charge to p and the ellipsis represents terms of order higher than 3 in {v/c). Notice that all 
functions are evaluated at present time. For a circular orbit, the term of (0) inside braces is 
of order (v/c) 2 times the Coulomb term and the term with the third derivative of position 
is of order (v/c) 3 times the Coulomb force. A detailed expansion, correct to fifth order in 
(v/c) is calculated in WM. The magnetic field caused by the charge q at p has the following 



series 

B = -|[/3xn] + ... (3) 

This first term is the Biot-Savart term and the next term in the series would produce a force 
of fourth order in (v/c) and along the normal, and it is not important for the present work. 
In this work we consider only the above terms of the electromagnetic interaction, and also 
consider the relativistic correction of Newton's law for the electronic motion up to second 
order in (v/c). 



Let us now discuss electrodynamics with retardation and self-interaction in the special 
case of the helium atom: We recall that the circular orbit is defined as one in which the two 
electrons are in the same circular orbit but 180 degrees out of phase @. Along such orbit, 
the total force acting on electron 2 can be calculated using (|2|) and the self-interaction of 
electron 2 to be 

F .|! ( x 1 + x 2 )-7§(l->s), (4, 

where R is the radius of the circular orbit. In (gj) we have also added the Coulomb attrac- 
tion of the a particle and used the approximation that the a particle is infinitely massive 
and resting at the origin. If the two electrons are in the same circular orbit but in phase 
opposition, the force along the velocity cancels out (first term on the right of equation |]), 
demonstrating that the circular orbit is a possible periodic solution of the electromagnetic 
equations up to third order. 

Notice the appearance of the dipole term in ([|), 

D = -e(x! + x 2 ). (5) 

The total far-field caused by the three particles depends linearly on the quantity D defined 
by equation (||), up to quadrupole terms. If this quantity is zero, the orbit is not radiating 
in dipole. The fifth-order terms of the Page series force and the fifth-order relativistic 
correction to the Lorentz-Dirac self-interaction introduce quadrupole effects. These effects 
would be important only in a much longer time scale, of order T/(v/c) 5 . In section VI we 
show that the circular orbits decay in a time of the order of T/(v/c) 3 . Therefore, the effect 
of quadrupole terms is small during the orbits lifetime. We start the next paragraph by 
discussing this importance of quadrupole terms in detail. 

Throughout this paper, it is very important to keep in mind the orders of magnitude 
relevant to atomic physics. For example, a typical value for (v/c) is (f/c) ~ 10~ 2 . Typical 
values for the width of the spectral lines is of the order of (v/c) 3 /T, which is the inverse 
of a time to perform about 10 6 turns along a typical orbit f26j. In the classical model of 
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the isolated hydrogen atom, because of dipole radiation losses, the energy loss during this 
linewidth-time produces dramatic changes in the frequency of rotation and therefore a band 
of dipole radiation is emitted, not a sharp line |^],[27|]. We explore this in another publication 
p7| , and we mention here this classical argument of Bohr |29[ just to stress that it does not 



apply for the circular orbits of a two electron atom. Along circular orbits, because there 
is no dipole radiation, a sharp line can be emitted for motion along the stable manifold 
28fl , provided the orbit is otherwise stable in an intermediate time scale. The quadrupole 



power radiated is of size (E /T)(v/c) 6 , where E Q is the coulombian energy of the orbit and 
T its period. This power times the linewidth time T/(v/c) 3 results in an energy loss of 
E (v/c) 3 , which is consistent with a sharp variation of the emitted frequency. The circular 
orbits should decay to the ground state or ionize in a time corresponding to the inverse of 
the linewidth, and for us here it is only important the fact that they can radiate a sharp 
frequency in the process. 

A word of caution should be said about the fact that a term in the Page series of 
order higher than two represents a singular perturbation. Therefore, when one such term is 
included in the equations of motion, there will be solutions to the dynamics that are not a 
perturbation of a mechanical coulombian orbit. We call these solutions "non-mechanical", 
as opposed to " quasi- mechanical" regular perturbations. In this work we investigate only 
quasi-mechanical regular perturbations of circular orbits. 

If one wants the most generic "stationary-state" , where the helium atom does not radiate 
in dipole, then one must have D = — e^+x^) = for all times. If we integrate the condition 
Xi + x 2 = twice in time, we get 

xi + x 2 = a + bt. (6) 

The constant b must be zero if the electrons are bound to the center of force at the origin. 
Inspection shows that a must also be zero, which is necessary for the existence of a solution 
to the equations of motion with the Coulomb force terms. Along orbits that satisfy ([]) with 
a = b = , the interaction between the electrons just renormalizes the charge at the center 
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of force. The most general "quasi-mechanical" stationary orbit possible is then an elliptical 
orbit, with circular orbits as a particular case. In this work we do not consider elliptical 
orbits, but they could be studied in a way analogous to the circular orbits. The other kind 
of possible stationary orbits are symmetric collinear motions of the two electrons. These are 
orbits of zero angular momentum and they are also singular orbits (Coulombian solutions 
with zero angular momentum and zero dipole would fall onto the nucleus). We discuss them 



briefly in section [VII . 



If one is not interested in non-mechanical orbits, it is convenient to truncate the Page 
series to third order and assume that the truncated system describes all the essential electro- 
dynamics and that the next terms only introduce a small stochasticity. For non-mechanical 
orbits, such as zero angular momentum orbits, the Page series is not convergent and it might 
be necessary to keep all orders, like in Eliezer's theorem p7, 30~l,|31|l . In the case the series 



is at least asymptotic, the successive orders become important in successively longer time 
scales. For example, about a periodic orbit of period T, the second order terms produce 
deviations from the Coulomb dynamics in a time of order T/(v/c) 2 . The third order terms 
take a time T/(v/c) 3 to influence the dynamics, and so on. In this work we consider the 
well defined dynamical system obtained by truncating the Page series interaction to third 
order, and we include the self-interaction to third order in (v/c) as well. 

III. COULOMBIAN STABILITY OF CIRCULAR ORBITS 

In this part we consider the non-relativistic Coulomb dynamics of a two-electron atom 
in the neighborhood of a circular orbit. The plain Coulomb stability problem has already 



been considered by many authors using other approaches, and we do it again to 

introduce our perturbation scheme. The Coulomb interaction is order zero of the Page series 
interaction (Q), and the scheme we develop allows for the inclusion of higher order terms of 
the interaction in an easy way. An alternative equivalent way to perform this calculation 
is to transform to a coordinate system rotating with the frequency of the circular orbit. In 
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this system the circular orbit is a fixed point of an autonomous vector field, the Jacobian 
matrix is independent of time, and the parametric problem is replaced by a linear eigenvalue 
problem. The disadvantage is that one has to transform all the terms of the Page series to 
the rotating coordinates. We get back to rotating coordinates later on. 

Newton's equations of motion for the Coulombian two-electron atom are 

z <? < ^ ^ ■ 

K la K 2a 
Ze 2 e 2 

M e Xi = -— g-(xi - Xa) - ^-(X 2 - Xi), 
K la K 12 

Ze 2 e 2 
M e x 2 = --^(x 2 -x a )-- r (x 1 -x 2 ). (7) 

where x Q , xi and x 2 are the position vectors respectively of the nucleus, electron 1 and 
electron 2, Ri 2 = |xi — x 2 |, R\ a = |xi — x Q | and R 2a = |x 2 — x a |. In the special case of 
helium, the nucleus is an a particle, but we will use the index a to label the coordinate of 
the nucleus in the generic case as well . The nucleus has an arbitrary charge of Ze, and the 
electrons have charge — e. The circular orbit periodic solution of (|7|) is 

x a 0, y a 0, z a 0, 
x 1 — Rcos(ut), yi — Rsax{ut) , Zi — 0, 
x 2 = —Rcos(ujt),y 2 = — Rsm(cjt), z 2 = 0. (8) 

According to (0), the frequency of the orbit is related to R by 

M e ul = {Z-\)^. (9) 

In this section, u and u are the same, but we will see later on that because of higher order 
corrections, io Q is only the first term of uj in powers of (v/c) 2 . To simplify the notation we 
define the quantity 

1 



8Z- 2 



which specifies a generic two-electron atom. 
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Linearizing (|7|) about the circular orbit @, we obtain a parametric linear differential 
equation with coefficients periodic in time and period T = tt/uj, a.s we show in appendix 
If the Floquet exponents are all nondegenerate, one can find a complete set of solutions of 
the form p3| 

<5x a = exp(2iufit) x" exp(2ma;t), 

n 

<5xi = exp(2iu;/xt) x* exp(2ma;t), 

n 

5x 2 = exp(2iufit) x^ exp(2mi<jt), (10) 

n 

where the Floquet exponent fi is a complex number defined in the first Brillouin zone, 
— 1/2 < Re(fi) < 1/2. From now on, an upper index should not be confused with an 
exponent, and takes values a, 1, and 2 to label the nucleus, electron 1, and electron 2, 
respectively. Notice that for the Floquet components we use an upper index, but to label 
coordinates as functions of time, as in (|^), we use a lower index (to distinguish it from the 
Floquet components). To bring the variational equations to normal form, we define the 
coordinates 

^ I K ■ K\ K ^ ( rp K I „•„ K \ 

~ ~\/2 ~ = iy/2 y-nli 

where again the upper index is not to be confused with an exponent, and takes the values a, 
1, and 2. Next we calculate Hill's secular determinant [p3| |, which reduces to the evaluation of 



a 6 x 6 determinant in this case of circular orbits. As a simplification, let us define n = n + /i 
to be the running variable in the summations of (|T0|). Notice that n defines a frequency of 
linear oscillation in units of 2u, as of (H^)- Last, to introduce the physical intuition in the 
problem and explore the symmetries, it is convenient to define the coordinates x r and x^ as 

X r = Xi + X2 — 2x a , 

X d =Xi-X 2 , (11) 

respectively, which introduces two new letters, r and d to appear as labels. In helium, 
{Z — 2), the coordinate <5x r is the dipole moment, and a nonzero amplitude for this variation 
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will always mean that the normal mode radiates in dipole. The dynamics of the 5x<f variation 
is decoupled from the 5x a and 5x r variations about circular orbits. Before we write the 
equations, we need still another definition 



1 



V n = ^C-l ~ 4^ + l" (12) 



It is convenient to think of the quantities and V£ as the x and y components of a two- 
dimensional vector K£. At this point one should take a look at Appendix |A], where we write 
the variation of the useful functional forms in terms of the vector Floquet components. Let 
us start by considering variations along the plane of the orbit. Using the results of Appendix 
[A], the planar variational equations of equation (|7|) can be written most simply in terms of 
the Floquet components of fllPf ) as 

-4n 2 xi-^(x£ + 3K£) = 0, 
-4n 2 x° + 4Z0(x; + 3K;)/^ = O, 

-4n 2 x^ - 4Z0(1 + -)(x; + 3K;) = 0, (13) 
Q 

where q is defined as the ratio of the mass of the nucleus to the electronic mass. For example, 
in the case of helium g « 7344 . For the Coulomb dynamics only, also the last equation of 
(|TB|) depends only on the radiation coordinate. It is good to have a scheme in mind to keep 
track of what we have done so far: equation fll3|) is the variation of (0) divided by M e u; 2 , 
which we write schematically as 

M5A(°) - 5F(°) _ 

In part IV we consider the variation of the relativistic corrections to the acceleration, which 
we call 6MAW, and the Page series second order force, which we call 5F( 2 \ The equations 
of motion will then be written schematically as 

MSA™ - 6FW SHAM - 5F^ 

M^l + M^l + • • • " °- 
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In this way, we keep adding higher order matrices and calculating Hill's secular determinant 
up to an order. Since we will be adding the matrices, it is necessary to keep some convention 
about the order in which the equations appear. The convention is that we always repeat 
what we did in (|T3|), that is, first the equation of motion for <5xd divided by M e u 2 , then the 
equation of motion for the 5x a divided by M a uj 2 and then the equation for <5x r divided by 
M e u 2 . 

We recall that n — n + //, and there is one equation for every value of n. In principle 
this would lead to an infinite Hill's secular determinant, which is the case for an elliptic 
periodic orbit. For the case of circular orbits, when we go to the £ and \ variables, we find 
that the £ n variables couple only to Xn+i and vice versa. This is obtained by taking linear 
combinations of the x and y components of the vectorial equations fll3|), a procedure which 
we call "unvectorizing" . We developed this convenient mnemonic method of sketching the 
calculation to avoid making mistakes in an otherwise lenghty algebra (even though we did 
it with the symbolic manipulator Maple). The resulting matrix equations are 



-3i 
2 



3/ 
2 



Si, 



X n +1 



0. 



(14) 



and 



— n 2 





Z(j)y 


3iZ(p 









—n\ 


—3iZ(p 


Z(j)y 




Xn+1 








-(n 2 + A) 


-3iA 




Sn 








3iA 


-« + A) 




Xn+1 



(15) 



where n+ = n + 1, y = (1/ q) and A = Z(f)(l + 2y). Finally, let us solve fli~4l) and 
the roots n. For ( |H| ) the roots are 



for 



-1,0,-1/2,-1/2. 



Notice that the secular problems of (0) or ( |T5D involve only n and n + = n + 1. Once n is 
calculated from the secular condition, we can always recover n and [i because any complex 
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number can be decomposed in a unique way as n = n+fi with n integer and |-Re(/i) | < (1/2). 
After we find n, this defines that the only two nonzero components are the nth and (n+ l)th 
in ( |I0D for the normal mode oscillation. An eigenvector calculation should then follow to 
determine the ratio of these two Floquet components. 

As regards the eight roots of equation (ITH) , four roots are 0, 0, —1, —1. From the general 



theory of linear ODE's, at a double root like n = 0, the general solution is a "quasi- 



polynomial" linear function of time p3fl . The doubly degenerate root n = is then respon- 
sible for the homogeneous translation solution with a constant velocity of the coulombian 
helium. Evaluating the determinant of ([15]) and equating it to zero, one finds that the other 
roots are given by 

1 

n = -- + e, 



where e is a solution of 



e 4 + (2A - £) e 2 + (± - 2A)(± + 4A) = 0, (16) 

which is a quartic equation with solutions 

e 2 = \ - A ± v/A(9A - 1). 

Inspection of the above formula shows that there is always a pair of roots with negative e 2 and 
a real pair of roots. The pair with negative e 2 describes an instability, which was first found 
by Nicholson and in this work we refer to it as the Coulombian radial instability P,|23|. 
This is an exponential growth in a time scale of the order of one cycle. The corresponding 
roots are given by 



n = - l - ± i^A - \ + ^/A(9A-1), (17) 

where i = a/— 1. One also finds that <5x r is not zero for the radially unstable mode, which 
implies that perturbations along this unstable mode radiate in dipole. We discuss some 
consequences of this in part V. The other two roots of (0), with positive e 2 are approximated 
by 
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n 



2 V 4 



A(9A 



(18) 



Notice that the quantity ( that appears in the paper of Poirier is equal to 8 times our A 
(A = Z(l + 2y)/(8Z — 2)) and that we are including the dynamics of the nucleus as well. 
To recover the results of Poirier one should put y = in our definition of A. 

Last, we consider the oscillations perpendicular to the plane of the orbit, which we call 
the z direction. In linear order in the oscillation amplitude, this oscillation is decoupled 
from the oscillations along the plane. Starting from equation (^), and linearizing about @ 
we obtain 



— + bz d = 0, 



5z n 



5z r 



8Z(f)y5z r = 0, 



+ 8A5z r = 0. 



(19) 



The solutions to this linear problem are of type 



&z K = c K exp(2iu\). 



(20) 



The first equation is a simple separate linear equation and the solution to it with q / 
requires A = ±1/2. The next two equations of (0) become 



-4A 2 -8Z(f)y 
-4A 2 + 8A 

and a nontrivial solution requires A = 0, and 



A = ±V2A. 



0. 



(21) 



(22) 



The linear z oscillation is decoupled from the planar oscillation , but of course it couples 
at higher orders in the oscillation amplitudes. 

As we mentioned in the beginning of this section, in a frame rotating with the frequency 
of the circular orbit, the circular orbit itself is a fixed point of an autonomous vector field. 
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This vector field describes the Coulomb dynamics in the rotating frame. One finds that the 
inclusion of the second order and radiative terms still yields an autonomous vector field for 
the dynamics in the rotating frame. This is the reason why we were able to simplify the 
parametric equation down to a 6 x 6 linear system. In this rotating system one also has 
to diagonalize a 6 x 6 matrix (two plane coordinates for each particle), and the exact same 
problem along the z direction. 



IV. INCLUSION OF SECOND-ORDER TERMS 

In this Section we consider the inclusion of the second-order terms of the Page series and 
the second-order relativistic corrections to the electronic dynamics. It is known that the 
Lienard-Wiechert interaction truncated to second order in (v/c) is described by the Darwin 



Lagrangian [13 



L Darwin ^ X „ Wjlij | ~\~ TTij|Xj| ) 

i 2 8cr 

1 1 

~9 E T^ 1 - ^3 & ' *i + (*< • n<i)(*i " ^v)])* ( 23 ) 

z ij 'ij zc 

where the indices take the values 1,2, and a, = |xj — x 3 -| and n^- is the unit vector in the 
direction of x» — x,. This Lagrangian depends only on the scalar product of the velocities 
and the distance between the particles, therefore it is invariant under a global rotation of all 



the particles coordinates. By Noether's theorem |20| , this generates an angular-momentum- 
like constant of the motion associated with the symmetry, which is equal to the angular 
momentum plus a small functional correction of order (v/c) 2 . Because the Lagrangian is 
time-independent, there is also an energy constant of the motion. This makes four indepen- 
dent constants of the motion. According to a recent result on non-integrability of generic 
systems of ODE's, for an extra analytic constant to appear, some extra resonance condition 
must be satisfied by the linear frequencies 0. 

Next we determine the second-order correction of all the Coulomb eigenvalues. Let us 
start by correcting the frequency of the orbit as given by @. A given circular orbit is 
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characterized by the velocity \f3\, and from this one can calculate all the other quantities of 
the orbit: angular momentum, frequency and radius. The Page series gives a correction in 
powers of \/3\, and the first correction is of order \f3\ 2 , as it should be for any relativistically 
invariant dynamics ||16|| . Adding up the second-order magnetic (|3]) and second-order electric 
(0) forces acting on the electron along the circular orbit we find the normal force 

w - e2n _ l&| 2 e 2 n 
8R 2 ' 



Together with the second order relativistic mass correction, this determines the second order 
correction to the frequency of @ to be 



A simple way to include this second-order frequency correction is to multiply n and n + 
by 1 + (d> — |J \(3\ 2 on the zero-order matrices of (|i4|) and (|I~5| ) before adding the other 
second-order matrices. At this point one should look up appendix B, where we evaluate the 
variation of the second-order terms of the Page series, as well as second-order relativistic 
corrections. 

Adding up the second-order electric field, second-order magnetic force, the second-order 
relatistic correction for the electronic masses and the correction to the frequency, we obtain 
a perturbation to the matrices ( |l~4"D and ([15]). We write the equation for the Floquet com- 
ponents, which separates in two parts just like (|T3|) splits into ([[4]) and (|15|). The resulting 
second-order correction to be added to ([14]) is 



Piin) iQ 2 (n) 
-iQ 2 (-n + ) P 2 (-n+) 



Xn+1 



(24) 



where /V^) 



2) n 2 + 8d)n + 



3© 



and n+ = (1+n). 



Q 2 (n) = [(40 + 2) (n 2 -n) + f 
Let us first consider the correction in powers of \j3\ of the roots of (fH]). This is done by 
adding up matrices ( |T4"D and fl2"4]), and taking the determinant of the result. The result, up 
to second order in \f3\ is 
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(25) 



16n(l + n)(l + 2n) 2 + |,3| 2 {(7 + 60) + 8 (1 - 50) n 
+8 (3 - 130) n 2 + 16 (1 - 40) (2n 3 + n 4 )} + ... 

In the neighborhood of the simple root n = of ([14])? (H) takes the form 

16n+ (7 + 60) \(3\ 2 = 0, 

which yields n = _ il+MML. Analogously, in the neighborhood of n — — 1, ( P5| ) takes the 
form 

-16(n+ 1) + (7 + 60) |/3| 2 = 0, 

which yields n = — 1 + iZ±MML_ Last, and crucial for this work, is the bifurcation of the 
n = — 1/2 double root. About n = —1/2, (f25|) becomes 

-4(1 + 2n) 2 + 6 (1 + 20) |/?| 2 = 0, 

with roots 



B= _> ± |^3(l+M. (26) 

All the above calculations can be done by hand, and we checked it using the program Maple, 
version 4.0. Notice that the correction of the degenerate root in comes with a linear 
power of \/3\, differently from the simple roots, that are corrected only at order |/3| 2 (exactly 
because of this degeneracy). This root undergoes the fastest change of all for small (v/c), 
and one should expect it to be the first to accommodate a "new" resonance condition. In 
the construction of a resonant normal form constant, this frequency allows the resonance 
condition to be satisfied with the lowest possible value of \/3\. (In the spirit of Section VI, 
one easily finds by inspection that resonances among frequencies corrected only at second 
order lead to relativistic orbits, not very interesting for atomic physics). 

Next we calculate the second-order corrections of (Tl5|): Adding up all the second-order 
corrections to (TT5|), second-order electric, second-order magnetic, second-order relivistic and 
correction of the frequency we find the following matrix to be added to (|i~5"D 
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(2 - 40) n 2 

(2 - 40) n^. 

i2(n) iS(n) 

—iS(—n + ) R(—n + ) 






T{n) iW{n) 
-iW(-n+) T{-n+) 



x 



ca 



Sn 
Xn+1 



(27) 



where R(n) = -8(l + 30)n 2 + 80n, S{n) = 4 (1 - 0) + (20 - 4) n, T{n) = 
—2 (1 + 80) n 2 + 40?7, , and Win) = 2 (1 — 20) n 2 — 2n, and we have left out terms pro- 



portional to y = 1/q. To calculate the corrections to the roots of (|i~5|), we add ( p7|) to (|15|), 
take the determinant and equate it to zero. Here it is necessary to evaluate the determinant 



using Maple, because we are dealing with a four by four determinant. Adding matrix (|27|) 
to ( |i~5[ ) and taking the determinant we obtain, up to second order in \(3\ 2 



n 



256rT (1 + n) 2 {(n 4 + 2n A + (1 + 20) n 2 + 2Z<pn + Z0-8Z 

+ i W l> [ «' + »» + g + f + •§*)-■ 



^2 



'l + £ + ^ fl 

v 10 5 3 J 



(28) 



Notice that according to equation ( |2"7| ) we still have the degenerate roots at n = and 
n = —1. An inspection in (|28|) shows that from all the roots given by (|i~5|) , only n = and 
n. = — 1 are still degenerate after the inclusion of the term in \/3\ 2 . We will not develop it 
here, but for y ^ these degenerate roots become simple roots and each other root and 
—1 is corrected by terms proportional to y\/3\ 2 . The roots of (|T8| ) are nondegenerate and are 
corrected only at second order as 

1 



n 



±\ A 



4 



+ V'A(9A-l)±c(Z)|/?|< 



(29) 
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The function c(Z) is a complicated function of Z, and we give its value for the most 
interesting values of Z, namely C(2)=0. 15241, C(3)=0. 11926 and C(4)=0.10126. 

As regards the second-order corrections for oscillations along the z direction, the equation 
for Szd is changed to 



§+^+im(~-20)§ 



0. 



and the roots are changed to 



A = ±i[l-(^-^)|/3| 2 ] 



In an analogous way, we find that the second order correction to (pip is 



4(±-0)A 2 



(30) 





16(2Z0-|-0) -120A 2 

Adding (^) to Q2TJ) , taking the determinant and equating it to zero, we find that the A = 
double root is preserved and the roots of 



are corrected to 



(31) 



This completes the calculation of the second-order correction to all the coulombian roots. 
In the next section we calculate the third-order corrections, and table 1 shows all the roots 
corrected to third order. 



V. INCLUSION OF THE DISSIPATIVE THIRD ORDER TERMS 

The terms of order higher than two in (|2|) are singular in the sense that they introduce 
the third derivative, and bring up a new solution to the dynamics. The next natural step 
in the study of the stability of the circular orbit would be to include the third-order terms 
in the calculation of the eigenvalues. We can still use resonant normal form theory in the 
presence of complex eigenvalues, and try to find a new resonance condition, which would 
then produce a constant of the motion up to radiative terms. We will find that many of 
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the resonance conditions that are satisfied by the Darwin interaction are destroyed when 
the third order terms are included. For example, resonance fl3"5|) is not satisfied anymore, as 
one can easily check using the imaginary parts as listed in table 1. This will not come as 
a surprise, and it signifies that the radiation makes those second-order stable orbits decay. 
If this is the cause of the decay, it is natural to expect that the time of decay should be 
of the order of T/(t>/c) 3 , which is actually the case in atomic physics. The introduction of 
the radiation is then seen to do two things: first, it makes the orbit to "dive" to the center 
of the resonance island in a slow time scale, emitting a sharp frequency. Second, since the 
resonance condition is not exactly satisfyed, the electrons can escape from the island during 
the long "dive", and decay to a lower energy state or ionize the atom. Another cause of 
the decay is the nonconvergence of the resonant normal form, whose optimal truncation 



produces a constant only for a finite time scale |34j . 

We now calculate the third-order corrections, in a way analogous to the second order. 
Starting with the third-order correction to ([H]), we find the third-order matrix 



.256i0|/3f 



A 3 ) 








c a 


(2-Z) 1_ 







which describes the radiative correction to be added to (|2~ID , in disregard of terms propor- 
tional to y\/3\ 2 . Solving the perturbed secular determinant, we find that the A = double 
root is preserved and the roots of (|22|) are corrected to 



A = ±4(2ZJ)[1 - h\p\ 2 ] 



+f*Vl/3| 3 , 



(32) 



where again i stands for the complex unit i = y/—l. 

Last, let us calculate the third-order corrections to (|14D and (|15|). Because of symmetry, 
one finds that there is no third order correction to ([14]). This is readily seen because the 
third-order variational force is the same for both electrons, being proportional to the <5x r 
variation. The third-order correction to (ITS]) can be calculated using the third order of (El), 
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in an analogous way we used for the second order. We find the following matrix describing 
the third-order correction to ( |T5| ) 



128 







2(2 - Z)n 3 2n 3 

2(2 - Z)n% 2n\ 



x 



Cot 
S n 



Xn+1 



Sri 



Xn+1 



Adding the above matrix to (15), add the second-order (27) and take the determinant. The 
result produces a correction to equation ( p8| ) given by 

16384|/3| 3 i 



5 3 = ^-(1 + nyn\l + 2n) x 

(n 4 + 2n 3 + n 2 + Z<pn 2 + Z<pn + Z<f>). 



(33) 



From the above, it is easy to see that the eigenvalues n = —1/2, n = — 1 and n = do 
not acquire any imaginary part at third order. It is also straightforward to add (|33|) as 
a perturbation to (pj|) and calculate the imaginary correction to (|18|) (using the program 
Maple). The imaginary correction is of type id(Z), and the most interesting values of d 
are d(2)=0.766, d(3)=0.740 and d(4)=0.731. In table 1 we show all the corrections of the 
eigenvalues, including order \/3\ 3 correction. We only showed in the table 1 the correction 
to the regular roots, by which we mean the ones that were already roots of the Coulomb 
dynamics. The introduction of the third order brings up nine more singular roots, which we 
do not consider. Last, we did not show in table 1 the unstable pair of eigenvalues (|H]) of the 
radial instability either, which are part of the eighteen coulombian eigenvalues but are not 
interesting for resonance conditions. The unstable mode should not take part of a resonant 
condition, but of course it will be part of the next orders of the normal form. 
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VI. SOME SPECIAL ORBITS RESULTING FROM THE RESONANCE 

CONDITION 

In this section we explore the orders of magnitude of some of the orbits. A rigorous 
approach to the material in this section would be to study the nonlinear stability first. This 
very hard task is not done yet and in appendix C we outline what remains to be done. The 
section is designed to inspect some of the stable orbits predicted by resonant normal form 
theory using the information already at hand. The guiding dynamical principle we use is 
that resonances with the minimal integer multipliers are the most important. It is known 
in general |35 that the size of the resonance islands varies as exp(— o), where o is the order 



of the resonance. This solves the paradox as to the infinite number of possible resonances: 
the ones with a high order occupy an exponentially small area of phase space, which makes 
them very unlikely. The situation is analogous to quantum mechanics, where there is always 
an infinity of energy levels. Nevertheless, in practice only a very small finite number of them 
exists in experimental situations. For example in the Balmer series, only the first twelve 
frequencies of the series can be observed as emission lines in very diluted gaseous states 
|19| , |2~9|1 . Since very high quantum states are too extended in space, one needs very rarefied 
gases and large astronomical masses of gas to produce a measurable signal. We consider in 
this section only some of the resonances with minimal ordering. The comparisom with the 
quantum results is very illuminating because it immediately suggests what further results 
are worthed pursuing with nonlinear dynamics techniques. It also brings up some of the 
strengths (surprisingly good agreement), and weaknesses (lack of a selection rule so far) of 
the research so far. In principle the dynamical system with retardation and self-interaction 
is a nonlinear system that can emit sharp lines, but it does not have to agree quantitatively 
with the quantum results. It is very intersting that it does so to some extent. We have by 
no means exhausted the material that would fit in this section. 

An approximation that we use in examining the resonance conditions is that we only in- 
clude the correction to the frequency that comes with linear order, because this is the largest 
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correction for small Let us now consider the special case of helium: As we mentioned 
before, according to the Darwin Lagrangian, helium has always an angular momentum-like 
constant of motion and an energy-like constant. In the neighborhood of some select circu- 
lar orbits, another complex analytic constant might exist, which could make those orbits 
nonlinearly stable. The constant we find here involves the amplitude of the normal mode 



corresponding to (22) in a combination with normal modes along the plane. This is possible 
because these modes are coupled at higher orders in the oscillation amplitude. We do not 
want to involve the linear modes describing the circular instability in the first resonant term 
of the constant, but the amplitudes of the unstable mode will naturally appear in the higher 
monomials of the series for the constant. By the way, this is how the circular instability 
can be equilibrated inside the small resonance island defined by the constant. Besides, the 
circular pair ( |I7|) are complex eigenvalues and would not satisfy a simple resonance condi- 
tion in combination with the other real eigenvalues, (one would have to vanish the real and 
imaginary parts of the resonance condition separately, which would require two integers at 
the best.) 

The circular orbit is a fixed point of the autonomous vector field describing the dynamics 
in a system rotating with the frequency of the circular orbit. To develop the normal form 
about the fixed point, we must move to this rotating frame. We recall that n is a frequency 
of oscillation, in units of 2u, of the variational dynamics, according to fllQl). In the rotating 
system, the new frequencies for oscillation along the plane are found by adding 1/2 to the 
formulas ([18]) and (p6|). The frequency ( f2~2|) for the z-oscillation is unchanged. A new 
resonance condition involving the root fl26|) is the easiest to be satisfied for small values of 



|/3 1 . Therefore, we suggest to look for a resonance among the frequencies of (0), fl2"2"|), and 
(p6|), which we rename as u>i, u>2, and u^, respectively, in the rotating frame, and which in 
the case of helium evaluate to 

U! = ~ 0.5345... 



/3 + V32 
C02 = y 2 g ^ 0-5560 
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12 

w 3 = yy|/3|^ 1.3093|/3|, (34) 

in units of 2u, and we have disregarded corrections proportional to y and \(3\ 2 . According 
to standard normal form theory, the necessary condition to have an additional analytic 
|||7| constant of the motion in the neighborhood of a fixed point is a resonance among the 
frequencies. By inspection, we find that a new quartic resonance involving the above three 
frequencies and with the minimal integer multipliers is of type 

oj\ — u>2 ~\~ 2u 3 = 0. (35) 

This resonance is satisfyed for |/3] given by 

\P\ = 0.0082.... 

The Coulombian binding energy of a circular orbit in a two-electron atom can be written as 
E = —mc 2 \{3\ 2 . This energy, in atomic units, for the above value of \(3\ is 

E = -1.265a.«. 

Of course, other resonance conditions are possible: for example, we could put an integer 
number in the resonance condition, as in 

Ml — U) 2 + 2t2Cl>3 = 0, (36) 

and the corresponding values of are given by 

n 

The discrete circular orbits corresponding to this have binding energies given by 

1.265 
E = —a.u. 

n z 

In the next section we show how to construct the extra complex constant of the motion in 
a Maclaurin series (resonant normal form) when condition ( p6|) is satisfied. 
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The surprising fact that we were able to pick particular circular orbits is a genuine 
signature of the nonlinear dynamics, because the linear eigenvalues depend on the circular 
orbit through the parameter \/3\. The extra constant of motion, together with the four other 
constants of helium can stabilize the orbit by returning the unstable direction back to the 
neighborhood of the special orbit. This will happen only inside some resonance region | |24| . 
We discuss this further in appendix C. The frequencies of the resonant orbits satisfying ( |36|) 
are given by 

0.8101 
oj = 5 — a.u., 

n A 

and the frequency of the z-oscillation in the stable manifold of the orbit can be obtained by 



multiplying the above frequency by 2y2/7, as of ( |22"D 

0.866 



w. 



n 3 



(37) 



Supposing that the dynamics in the neighborhood of the resonant orbits is stable , one 
could expect that the stable oscillations about this orbit could emit a sharp line. The only 
condition to emit a sharp line is to oscillate with the same frequency for a long enough time 
(of the order of the inverse of the width of the line). This condition is fulfilled because of the 
finite time stability of the resonant orbit. The correct frequency of these stable oscillations 
can only be obtained after the nonlinear stability is performed, by linearizing about the 



center of the resonance islands. Here we will assume that the frequencies of |37] are an 
approximation to those. 

For the circular orbit corresponding to (|37|) with n = 1, the frequency of the z-oscillation 
in the stable manifold is 0.7956 atomic units. The transition from the first excited state 
of parahelium to the ground state (2 X P — > l 1 ^) corresponds to a frequency of 2.9037 — 
2.1237 = 0.7799 atomic units [PT] , which is a 9% difference. For n = 2, ( p7|) evaluates to 
0.1083 atomic units, and the frequency for the transition (3 1 P — > 2 X S) in parahelium is 
2.1459 — 2.0551 = 0.0908 atomic units, which is again a 9% difference. Last, the asymptotic 



form of the quantum energy levels of helium [18| , both parahelium or orthohelium, is 
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E _-z 2 {Z-lf 

2 2(n r + L + l) 2 ' 



with Z — 2, which is a first approximation to the Rydberg-Ritz spectroscopic term []T9 



The frequency of the line emitted by transitioning to the neighboring level, calculated from 
the above formula with AL = 1 is approximated by 

1 

w = 



(n r + L + l) 3 ' 

This above equation is a quantum formula which we write just to compare with 
Notice that (^) agrees with it to within 13%. 

We can also produce an estimate for the width of the line as follows: The third-order 
correction to the frequency of ( p2|) (w z ) is imaginary and with the same sign for the two values 
of A. Therefore the ^-oscillation is stable and decays with a coefficient which, according to 
([2C)|), is the imaginary part of 2a; A. For example in the case of helium, Z = 2, and 4> = 1/14 
, this imaginary correction to (|37| ) evaluates to 

and this oscillation is then part of the stable manifold of the circular orbit. Along this 
decaying oscillation, the perturbed circular orbit decays back to the perfectly circular spe- 
cial circular orbit. The radiative self- interaction has already been shown to produce good 
approximations for the linewidth in other situations ||26|| . If we evaluate the above result 
for the line at ui = 0. 7959a. u., we find that it is 1.5 times the experimental linewidth for 
this transition []T7] , |3"5(] . We stress again that this is not exactly the probability of decay of 
the orbit: As we already mentioned, the resonance condition Q3"5] ) is not satisfyed with the 
inclusion of the third-order terms, which implies that the constant of motion is destroyed. 
Because of this, the perturbed orbit can decay not only back to the perfectly circular circu- 
lar orbit, but also to the lower-energy ground state. The dynamical linewidth would be the 
full probability to escape from the atracting resonance region of the stable orbit, and the 
calculation of this is beyond the scope of the present paper, but since the damping of the 
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oscillation is causing the decay, one would expect a number of the order of this damping for 
the inverse of the linewidth, which is again a good agreement, since we found 1.5 times the 
correct quantum result. 

Let us briefly consider the case of the Li + ion: The frequencies of can be easily 



recalculated for the case of Li + by using the material of sections III and IV with Z = 3 and 
= 1/22 .We find 

out = - 0.5222... 



w 3 = J^\P\ - 1.2792|/3|, (38) 
and a simple resonance condition like 

U)\ - tu 2 + lo 3 = 0. (39) 
will determine the value of \/3\ to be 

\/3\ = 0.0125. 
The stable orbit has a Coulombian energy of 

E = -2.932a.u. 
and the frequency of the z-oscillation mode is 

w z = 1.97a.u. 



again in very good agreement with quantum mechanics. The first two quantum energies of 
para-lithium Li + |3J| are: ground state, l 1 ^ : E = — 7.278a. u. and 2 1 P : E = — 5.30a. u. The 



frequency of the dipole transition becomes w = 7.278 — 5.3 = 1.978a.u., in very good 
agreement with our above value of w z (within one percent). Notine that this time we 
started the resonance condition with n=l instead of n=2 as for helium. It would be futile 
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to go on without knowing the reason for the special integer combination. In appendix C 
we show that this selection rule requires the knowledge of the next term of the resonant 
normal form . This would be also the case if we tryed to apply this to find some states of 
H~(Z = 1) : It is known that H- minus has very few, if any excited states, and accordingly, 
the nonlinear stability should predict that all most of the circular orbits are unstable. The 
frequencies we obtained are in good agreement with the sharp lines of helium. The values 
of the energies do not agree so well with the quantum energies, but are nevertheless in the 
same order of magnitude. This comparisom should be repeated after one knows more about 
this dynamical system. 



VII. DISCUSSIONS AND CONCLUSION 

Historically (1912), coulombian many-electron atoms ("saturnian atoms") were investi- 
gated by astronomers |32j prior to Bohr, who was well aware of these studies. For oscillations 
perpendicular to the plane of the orbit, the Coulomb dynamics is stable and the ratio of 
many lines obtained by Nicholson for perpendicular oscillations agreed with the spectra of 
the Orion nebula and the Solar corona f32|. Of course Nicholson was assuming a special 
radius for the orbits, which he did not know how to calculate, and this radius would disap- 
pear when one took the ratio of two lines of the stable manifold of an orbit (because of a 
degeneracy of the Coulomb interaction). Bohr was originally favorable to the use of ordinary 
mechanics to describe those stable perpendicular oscillations |29|] . For oscillations along the 



plane of the orbit, Nicholson first found the now well-studied Coulombian radial instability 
PJ22||, which was then a hindrance for the theory |29j. It is of historical interest to stress 



that it was the radial instability that first motivated Bohr to postulate a discrete set of 
special, more stable orbits, which proved to be a very fruitful intuition [^5]. The original 



critic of Bohr |29[| to Nicholson ||32j| , was that the circular instability would make the atoms 
too "fragile" to disintegration, and unable to emit a sharp line. This led Bohr to conjecture 
that, along some special stable orbits, some "quantum" mechanism would supersede normal 
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mechanics and prevent the radial instability |29| . 



As we mentioned in the introduction, the classical chaotic Coulomb dynamics of helium 
has been widely studied in connection with the recent interest in periodic-orbit quantization 



22|,|23"|j. The circular instability is now known to lead to self-ionization, after a long-term 



chaotic transient, for most initial conditions 12311. Of course, those studies are for a different 



dynamical system, with lesser constants of motion, once the extra constant needs the retar- 
dation to exist. It would be of interest to repeat these numerical studies with the inclusion 
of retardation. 

As regards to where the circular orbits ultimately decay to, we conjecture that the lowest 
energy bound-state in helium could be a "non-mechanical" orbit of symmetric collinear 
motion. This orbit has zero angular momentum and zero electric dipole moment. The 
avoided three body collision can be provided by a singular mechanism analogous to the 
one of Eliezer's theorem for hydrogen |]30],[3TJ. In hydrogen, Eliezer's theorem predicts that 
the electron will runaway from the atom, which is closely related to the fact that one has 
always dipole radiation in atomic hydrogen. The symmetric collinear motion in helium is 
not radiating in dipole, and therefore one could expect a physical solution, differently from 
the case of hydrogen |3TJ . It is intersting to compare these zero-angular- momentum singular 



solutions of Eliezer's theorem to the divergent series of the Lamb shift |l8|Jl9|] in quantum 



mechanics. The Lamb shift appears in quantum mechanics because of a singular interaction 
with the electromagnetic field and is only pronounced for states of zero angular momentum 
fIT||T9|| . In the dynamical approach, the solutions with zero angular momentum and zero 
dipole will be non-mechanical collinear orbits that get very close to the nucleus. This in 
turn causes the Page series to diverge or be asymptotic at the best. 

A recent use of resonances in perturbation theory worth mentioning here was on the 
problem of the time of stability of integrable tori. Here one exploits the resonances among the 
unperturbed frequencies . Those results go under the name of Nekhoroshev bounds 

for Arnold Diffusion. Simply stated, the results say that the actions of an e-perturbed 
Hamiltonian system are kept approximately constant for a time of the order of 
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T~(±)exp ([!]<*), 

where a = (1/ef) and c? is the maximal number of unperturbed frequencies linearly indepen- 
dent over the rationals [f40| , f45|j . Every time there is a resonance among the frequencies, d 



is reduced and the torus has an exponentially longer time of stability. This phenomenon 
is named "stability by resonance" [fS|. In connection with these modern Arnold Diffusion 
results, it is interesting to mention that they shone new light onto the old problem of the 



ultra-violet catastrophe, which was the historical motivation for Plank's hypothesis |f40| , [46 
Today, also many numerical results exist showing that a set of coupled oscillators might 
never reach equipartition f47f , the reason being a super-slow Arnold Diffusion. 



As regards further research to be done, we have only been able to study the circular 
stationary orbits, and it would be of much interest to study the most general elliptical sta- 
tionary orbits and associated energies and spectra, which should present a richer structure. 
Along elliptical orbits one can still use regular perturbation theory, in the same way used 
for circular orbits here, but now the parametric problem associated is more complicated. It 



might be that the use of the Kustanheimo transformation |48[ to regularized coordinates 
will simplify the problem of elliptical orbits. For zero-angular-momentum collinear orbits, 
it might be necessary to introduce the retardation effects in a non-perturbative way. To 
"quantize" those zero-angular-momentum orbits using resonant normal form theory is an 
open problem of much theoretical interest. 

One would expect some discussion about spin: Notice that our simplified dynamical 
system is based on classical pointlike charges with no spin. Of course, we could include spin in 
the dynamics in a phenomenological way, similar to the usual way it is introduced in quantum 
mechanics. We have not done that yet. Second, quantum mechanical spin is a relativistic 
effect and introduces a correction comparable to the second-order retardation correction 
rSyiPfl . As a matter of fact, the Breit operator for helium is actually produced starting from 



the Darwin Lagrangian |18|,19| . We believe that a correct discussion of spin issues can only 



be made after we know more details about our dynamics, such as the frequency correction 
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due to motion along resonance islands and etc. 

As a summary, we presented a complete account of the linear stability of a two-electron 
atom along circular orbits in the presence of retardation and self-interaction. We calculated 
all the linear eigenvalues up to third order in (v/c). We considered the necessary condition 
for an extra constant and showed how to construct this extra analytic constant. We have also 
shown that electrodynamics assigns a dynamical system to helium, which can be responsible 
for the emission of sharp spectral lines. (The sharp line is emitted by "diving" to the center 
of the resonance island, radiating out some energy and decaying to a lower state). We found 
that the frequency of the sharp line associated with the lowest stable orbits agrees with the 
highest frequency in the spectrum of helium to within 8%. We do not know of prior results 
on the existence of these stable electromagnetic orbits, which appear as a genuine effect 
of the nonlinear dynamics prescribed by Maxwell's electrodynamics to atomic physics. We 
have barely touched the study of this dynamical system, and much research remains to be 
done, specially in the nonlinear stability along the lines of appendix C. 

APPENDIX A: VARIATION OF THE COULOMB INTERACTION 

In this Appendix we calculate the variation of the terms appearing in the Page series, 
equation @. To obtain the variation of the Coulomb force along the plane of the orbit, one 
has to evaluate 



where r = |x| and (Sx, 8y) are the functions of time representing the variation about the 
circular orbit 




5x 3(x8x + ySy)yi 



(Al) 



x = ±rcos(o;t), y = ±rsin(W). 



(A2) 



Substituting (^) into ( |A1| ) we obtain 

X 3 

5 (^)n = {8x - -(1 + cos{2ut))5x - 



sm(2ut)5y} 
sin(2ujt)5x} 



8{^) n = {8y-^{l-cos{2ut))8y-^ 



(A3) 
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where the subscript n indicates the nth Floquet component. Using equation (|T0D we find 

(5x) n = x n and (Sy) n = y n and 

x 1 

^)n = -^- 3 (x n + 3V n ), (A4) 

where we have used U n and V n as defined by equation (|12"1). An economic way to write this 
equation is 

^)„ = -^(x n + 3K n ). 

When we consider the second-order terms of the Schott-Page force and the relativistic cor- 
rection of the mass, we need to evaluate the variation of many other functions besides x/r 3 . 
This is done in Appendix |B[ in a way analogous to the above calculations. Last, along the 
z direction the variation is simply 

Xf z \ 5z 



APPENDIX B: VARIATION OF THE SECOND-ORDER FORCES 

In this Appendix we consider the terms of order (v/c) 2 in the force. One contribution 
comes from the first relativistic correction to the electronic masses (the a particle is at rest). 
The x component of this term is 

5M e aW = ^5(x ■ x)x e + ^5(|x e | 2 x e ) 

Substituting the equation for a circular orbit, we find 

M e u 2 ,„, 0r/ „ ,„ x Jx 



W e a?) = - cos(2o;i))^ 

z, to 

- sin(2^)^% + 2(sin(2cut) — - cos(2cut)^)l 

CO 2 CO CO 

Using the above equation, equation ( |TU[ ) for <5x and definition ( |I2"D for U n and L n we obtain 
the Floquet component 
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(SM e a {2) ) n = M e u 2 \(3 e \ 2 (-4n 2 x e n + 2n 2 U e n + 2inV e ) 

The y component can be calculated in an analogous way. To write the two components in 
a concise vector form we define an extra useful vector quantity 



V„ 



This quantity is obtained from the vector K (defined just bellow equation ( |T2| ) by exchanging 
the components and changing the sign of the second component). Making use of the above, 
we write the relativistic second-order force as 



(5M e a (2) )„ = M e u 2 \(3 e \ 2 (-4n 2 x. e n + 2n 2 K e n + 2mK e ) 

Using a procedure analogous to the above, we can obtain the variation of the second-order 
electric interaction of (pj). First, the variation of the second-order electric field produced by 
the electron at the point p on the observation point o is calculated. The vector pointing 
from the electron to the point p (where another particle is) is x = x& — x e (k designates 
the other particle). In all cases of interest we have x = — gx e . It is easy to verify that for 
electron-proton interaction g = 1 and for electron-electron interaction g = 2. According to 
(0), the variation of the electric field is proportional to the product of the charges times 

^ |&| 2 x $ e (n-/3 e )n 
2r 3 2rc 2rc 2 

^{(^Xxn + 3K n ) + 2nV (x^ + K£ ) + 



2ign{± e n ~ K) h 



(Bl) 



where r is the interparticle distance and R the radius of the electronic orbits. The tilde over 
x means the vector obtained from x by exchanging the components and changing the sign 
of the second component. 

Vn 
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where again the superscript k can take the values l,2,a,r and d, to represent a particle's 
coordinate or a combination of the coordinates, as defined in (|TT|). The second-order electric 
field produced by the nucleus on each electron is 

the same for the two electrons. Last, we calculate the variation of the second-order magnetic 
forces. The a particle is at rest and does not produce any magnetic field. The two electrons 
produce a magnetic field at the center of the orbit, and the Lorentz magnetic force over the 
nucleus is 

_ 4Zie 2 ul z 

where l z is the ^-component of the angular momentum of one electron. The variation of the 
electron-electron magnetic force done by electron 2 over electron 1 is 

e 2 \B I 2 

(*F2i)„ = l^L{2<n(2x^ - if 

+Kf)-(x d n + K d n )}. (B2) 

For the force done by electron 1 over electron 2 just interchange the indices and the sign in 
front of and K„. 

APPENDIX C: CONSTRUCTION OF THE EXTRA COMPLEX CONSTANT OF 

MOTION 

To construct the resonant normal form, it is again convenient to move to the rotating 
system. For this part we disregard the acceleration of the nucleus, so that x a = 0. It is 
also convenient to use the Hamiltonian formalism associated with the Darwin Lagrangian, 
written in terms of the coordinates x r and x rf , as defined by equation (|TTJ), and the conjugate 
momenta p r and p^. We make a last transformation to the normal mode coordinates of 
the Jacobian matrix about the fixed point corresponding to the circular orbit. Let U{, 
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i = 1, . . . , 12, be the coordinates corresponding to the normal modes of frequencies Ui, 
i = 1, ... ,12, respectively. These coordinates will be complex- valued linear functions of the 
coordinates and momenta. Because of the simplectic symmetry, the eigenvalues exist in pairs 
lu and — uj [^5|. We order the modes such that uyi-i is the coordinate corresponding to the 
frequency = —uJi- Because there are six coordinates and six conjugate momenta, we 

have twelve normal coordinates. For example, the normal mode solution for the coordinate 
Mi will be given by 

ui(t) = «i(0) exp(iuit), 

and so on for the other coordinates. We define these coordinates to describe the linearization 
about the fixed point corresponding to the circular orbit, so that at the circular orbit all 
coordinates are zero. 

We want to find an extra constant with an analytic form, which implies it has a formal 
Taylor series about the fixed point (not necessarily convergent). For instance, the integer 
powers appearing in this Taylor series are the cause of the integer numbers in the necessary 
resonance condition J7|. Before we go on, let us introduce some definitions to simplify the 
exposition: Because we consider a Maclaurin expansion of the constant of motion, some 
notation concerning monomials is in place: We write a monomial in the coordinates Uj, 
j = 1, ... 12 as 

where k is a twelve-dimensional vector of integer components. For example, if the ith 
coordinate is not present in the monomial, then ki = 0. We use the convention that when ki 
is negative, the coordinate to be used is the complement of ki to 12 with the positive power 
(fci2-j = — ki). For example 

m(t)- 1 EE U 7 (t), 

and so on for all the other coordinates. Mode 7, by definition, has a frequency that is the 
negative of u>i. With this convection, even though we are working with negative integers, 
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all the powers of the coordinates in the series for the constant are positive, as it should for 
any good Taylor series. 

Let us define the order of a monomial by 

o = \ki\ + \k 2 \ + • • . + \k 12 \, 

and the frequency associated with a monomial as 

At<; k = k 1 uj 1 + k 2 uj 2 + . • • + h 2 uj 12 . 

The necessary condition for the existence of the extra analytic constant is that some 
resonance condition must be satisfied among the linear eigenvalues. In the following we 
show how to construct the constant starting from a leading resonance 0-0. Let us assume 
that a resonance defined by k Q exists 

A^ ko = 0. (CI) 
Our proposed analytic constant must have a Maclaurin series expansion given by 

k 

where the Ck are constant numbers. This is an analytic function by construction, because 
it involves only integer powers of the coordinates. We now construct a constant where the 
lowest-order monomial is defined by k Q as of ( |C1|) , plus higher-order monomials necessary 
to vanish the time derivative at higher orders. For example, if the resonance condition fl35|) 
is our leading resonance, then we can construct a nontrivial analytic constant by starting 
the Maclaurin series with k\ = 1, kg = 1, k% = 2, and k 2 = k% = . . . = k\ 2 = 0, as of (|35|). 
To produce the successive powers of the constant, let us write the constant as 

C = U ko +^C kl f/ kl + ... (C2) 

ki 

where we have separated out the contribution of the leading term at k , the contribution of 
the monomials of immediately next order, which we index with ki, and the dots represent 
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monomials of higher than immediately next orders. This is a complex valued function, which 
means two real functions of phase space. Of course, starting the series with the lattice vector 
— ko will produce the complex conjugate of the same function. Next we evaluate the time 
derivative of flC2|). 

The time derivative of the leading monomial produces the same monomial multiplied by 
Aa>k = 0, plus higher-order terms, as 

-^(f/ k <>) = + £iV kl f/ kl + . . . 

at ki 

This monomial is then almost a constant, up to higher order monomials. To produce the 

higher monomials of the constant, let us focus on one such term of immediately next order, 

corresponding to a vector ki and with coefficient N kl . The time derivative of C is 

dC 

— = Au k0 C k0 U k0 + £ N kl U^ + ... 

ki 

+ ]TAcu kl C kl £/ kl + ... (C3) 

ki 

It is easy to see that the next order nonlinear contribution of the leading monomial to the 
time derivative can be canceled by choosing 

C kl = -iV kl /Acj kl . 

It is well known that this perturbation scheme does not produce a convergent constant, 
which is not a problem if we are investigating stability for a finite time scale. The usual 
procedure is to terminate the series using some optimal truncation p4 |. This produces a 



quasi- const ant for a long time-scale, which still provides stability in the finite time for the 
special orbits. The detailed investigation of this issue is beyond the scope of the present 
paper. Here we are mainly concerned with exploring the necessary condition for the existence 
of such a constant. Just for illustration, let us assume that optimal truncation means 
considering the first term of the series plus the largest next order term. As this term contains 
resonances, the largest next order term will come divided by the smallest resonance, which in 
our case is proportional to \/3\ 2 - To see what kind of constraint this constant can provide, let 
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us assume that the perturbations about the orbit have a magnitude of size (. The constant 
expressed in terms of £, with the first and most resonant next term included is 

A4+2n 

c-C p |/?|2 +••• 

In the above equation, the first power is because the order of the resonance in assumed 
to be 2 + 2n, and we are also assuming that the next most important term is of order 4 + 2n. 
The constant p will be a function of n and Z, and the selection rule will be p (n, Z) > O.The 
maximum excursion away from the unstable fixed point will be given by C = 

(2 



p(n,Z)' 

We see then that the sharp resonances can make this excursion small, which is consistent 
with a sharp frequency emited by the tangent dynamics. 

Let us briefly discuss the shape of the level curves of such a constant: Notice that we 
have involved the frequency of the z coordinate oscillation mode in the resonances of section 
VI, which makes the level curves of the constant involve displacements along the z-direction. 
These level curves will then have the shape of a torus encircling the circular orbit, and with 
an excursion amplitude given by the above formula. When p(n, Z) > 0, the six constants 
intersect to confine a bounded manifold inside the resonance islands. In quantum mechanics 
this information is provided by the normalizability of the corresponding wave function, and 
the boundedness of the region inside the resonance islands seems to be the analogous concept 
here. We study the detailed stability problem in another publication [|36| 

I acknowledge discussions with R. Napolitano, L.Galgani , A. Carati, A. Giorgilli and D. 
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FIGURES 

FIG. 1. Table 1: Corrections to the sixteen stable regular roots, (f> = 1/(8Z — 2) and 
A = Z<f>(l + 2y). 



45 



